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1. Introduction 

As more computing power and better algorithms have become available over the years the 
precision in lattice QCD results becomes more and more relevant. In order to have small errors 
on lattice QCD results it is important to have good control over the systematics involved in lattice 
simulations. Besides large volumes, quark masses close to the physical point, and large statistics, 
the continuum limit plays a prominent role because of the so-called topological freezing prob¬ 
lem [1-3]. As the lattice spacing is decreased and, hence, the continuum is approached it becomes 
more and more difficult for conventional simulation techniques to tunnel between topologically 
distinct gauge sectors during the simulation. As a result, one faces large autocorrelation times for 
simulations at lattice spacings smaller than about 0.05 fm [4]. With the conventional approach 
large statistics and/or large spatial volumes would be necessary to obtain meaningful results. A 
quite different and new approach to reduce the large autocorrelation times for quantities which 
couple to those topological modes is based on the introduction of open boundary conditions [5] 
which are, in practice, imposed in the time direction. Topological objects can then flow in and out 
through the temporal boundary, and, in fact, it has turned out that this approach is computationally 
advantageous [6] with respect to an - otherwise necessary - increase in volume and/or statistics. 

Therefore, the CLS collaboration [7] has started a major effort to generate gauge field con¬ 
figurations with open boundary conditions aiming at large statistics and volumes as well as pion 
masses close to the physical point in order to get good control over all systematics. There has 
been generated a relatively large set of ensembles quite recently using Wilson fermions with 2-\-\ 
flavors [8]. We report on first measurements on some of these ensembles, in particular we focus on 
first results of the SU(3) octet baryon spectrum and also briefly address the scale setting, see Sec. 3 
and Sec. 4, respectively. Before, we will start with a short presentation of some details about the 
simulations. 

2. Simulation Details 

A non-perturbatively ^ (a)-improved Wilson action with tree-level Symanzik improved gauge 
action and 2 -|- 1 flavors has been employed. We have simulated at pion masses ranging from 
about 420 MeV down to about 130 MeV following the strategy adopted from the QCDSF collab¬ 
oration [9] where the sum of bare quark masses is kept constant along the chiral trajectory. This 
strategy not only simplifies the tuning process of the simulations but also has the advantage that 
one can extrapolate, e.g., the combination rather precisely to the physical point as it 

turns out that in leading order chiral perturbation theory this combination is proportional to the sum 
of quark masses, -|-/m^/2 oc TrM = and, in addition, the corrections to this relation stay 


run id 

HlOl 

H102 

H105 

ClOl 

DlOO 

mjt K. 

420 MeV 

350 MeV 

280 MeV 

220 MeV 

130 MeV 

mx « 

420 MeV 

440 MeV 

460 MeV 

470 MeV 

480 MeV 


Table 1: CLS ensembles at j3 = 3.4. Note that the DlOO ensemble has currently only low statistics and is 
not considered in the analysis presented in this contribution. 
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Figure 1: Effective mass plots of the nucleon. The horizontal length of the red bars in the plots indicate 
the fit range which was chosen to be [10,18], the vertical position indicates the fitted mean value and the 
width corresponds to the fit error. Note that for the ClOl ensemble « 220 MeV) the statistics is smaller 
compared to the other ensembles. 


small down to the physical point. As the scale to is known to depend only weakly on the quark mass 
(see also the upper left plot in Fig. 3) we use the dimensionless combination (/)4 = 8to(m^ + m^/2) 
and match our chiral trajectories with different lattice spacings at a single point along the symmet¬ 
ric line, i.e. we have fixed TrM, respectively 1 / tc,-, for each lattice spacing to match our target 

estimate (p 4 \m„j=m,, = 1-15 for niud = Note that, at the physical point, we obtain 

= 1.117(38) with = 134.8(3) MeV and mk = 494.2(4) MeV from Ref. [10], and 
\/8to = 0.4144(59)(37) fm from Ref. [11]. 

Currently, we have generated ensembles at three different lattice gauge couplings with j3 = 
6/g^ = 3.4,3.55,3.7. In the following we present preliminary data only for our ensembles at 
j3 = 3.4 which corresponds to a lattice spacing a « 0.086 fm. Note that also for j3 = 3.55,3.7 
other quantities have been presented at this conference [12, 13]. We have summarized the pion 
and kaon masses of the ensembles with j8 = 3.4 in Tab. 1. Note that the DlOO ensemble has 
currently only low statistics and is therefore not be taken into account in the analysis below. The 
lattice geometry for the ensembles HlOl, H102, and H105 is 32^ x 96 and for the ClOl ensemble 
48^ X 96 where the temporal extent of the lattices is typically larger compared to simulations with 
periodic boundary conditions due to the boundary effects specific fo open boundary condifions. All 
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Figure 2: Fits to the octet baryon masses of our preliminary data (right) compared to fits from Ref [9] based 
on data from Ref. [20] (left plot was taken from Ref [9]). 

ensembles have large spatial volumes ; for all ensembles trij^L > 4 is satisfied. The HMC 

trajectory has length T = 2 and for the ensembles under consideration more than 4000 Molecular 
Dynamic Units (MDUs) for each ensemble have been generated. 

For the generation of gauge field configurations the openQCD package [14] has been used 
which has several algorithmic improvements built in. Besides the Hasenbusch trick [15], improved 
integrators [16], a multi-level integration scheme [17], and a deflated solver [18, 19], the twisted 
mass reweighting [6] should be pointed out because it is the first time that this technique has been 
applied within a larger set of simulations. In the light fermion part of the action a twisted mass 
term is introduced in order to push the eigenvalues of the Dirac operator away from zero and, 
hence, increase the stability of the HMC simulation. 

To correct for the additional twisted mass term in the light quark action one has to reweight 
the observables. In addition, one also applies reweighting with respect to the strange quark action 
to account for errors coming from the rational approximation which is used to simulate the strange 
quark. The ensemble average of an observable {O) = then includes the reweighting factor 
W = WqWi with a light and strange reweighting factor Wb and Wi, respectively. We found that 
these reweighting methods work quite nicely in practice, more details can be found in Ref. [8]. 

3. Baryon Spectrum 

In this section we present first preliminary results on the baryon spectrum from simulations 
with open boundary conditions. In Sec. 3.1 we give some details about the computation of the spec¬ 
trum, in Sec. 3.2 we discuss the chiral extrapolation of our results and also compare our hndings to 
results from the literature. 

3.1 Computational Details 

For the computation of the baryon spectrum we have used standard relativistic interpolators, 
e.g. for the nucleon we have used In = EabcUa {ubCy^dc) ■ The source position in temporal direction 
was always chosen to be at the center of the lattice, i.e. t^rc = 47 where t S [0,... ,95], the spatial 
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source position was chosen randomly. We have been investigating smeared-smeared point-to-all 
correlators where we have used 100 steps of Wuppertal smearing on APE smeared gauge links 
for both the source and the sink. The fit range was chosen to be t G [10,18]. For the ensembles 
HlOl, H102, and H105 all available 4000 MDUs for each ensemble have been analyzed, for the 
ClOl ensemble 1000 MDUs have been considered so far. For the analysis a custom version of the 
CHROMA software package [21] has been developed where also the deflated solver [18] is used. 

In Fig. 1 we show the effective mass of the nucleon for the ensembles under consideration. For 
the three ensembles where all configurations have been analyzed we observe a nice plateau while 
for the ClOl ensemble one would expect that a less noisy plateau will develop once full statistics 
have become available. For £ and S the corresponding plots look similar. Note that it has turned out 
that for this preliminary analysis the smearing was not chosen optimally. We are currently running 
measurements with full statistics and optimized smearing and, hence, we expect to improve on the 
errors present in this study. 

3.2 Chiral Extrapolation 

In order to chirally extrapolate to the physical point we make use of SU(3) chiral perturbation 
theory. The octet baryon masses in terms of the pion and kaon masses to order ^{p^) are given 
by [22], 

niN = mo — AboM^ + Abf {m\ — — 2bo {2M\ + + • • • , 

m\ = niQ + —bi) (—4M|- + M^) — 2bo {2M\ H-, /-o i 'v 

j fj.tj 

= triQ — Abj)M^ — 2bQ (2M^ -\-M^ + • • • , 

ms = mo - AboMl - 4bp {mI -MI)- 2bo {2MI + ) + • • • . 

The average nucleon mass and average pion mass Xji are defined by X^ = {niN + + ms) /3 

and X^ = {2rn\ + m^) /3. In leading order chiral perturbation theory one finds, 

Xn = mo - 2bo (2M| + M^) +... . (3.2) 

Since 2m|-+m^ is approximately constant along our chiral trajectory also the average nucleon mass 

,„2 _^2 

stays constant to leading order along that line. We define the dimensionless quantity v = ^ and 

plot the fits of our data to Eq. (3.1) in the right plot of Fig. 2. The errors on the masses come from 
uncorrelated fits where we have performed a jackknife analysis with binning. The left plot of Fig. 2 
is taken from [20] where chiral fits have been performed to data from QCDSF [9]. Comparing 
those to our preliminary results we find nice agreement. 

4. Scale Setting 

To set the scale at a certain lattice coupling one would choose a dimensionful, (ideally) phys¬ 
ical quantity and compare it to the lattice result. In order to do so we have to, in our case, chirally 
extrapolate to the physical point. For the following discussion we define as the physical point the 
value of Ki = at which the ratio of strange to light ^(a)-improved AWI quark masses takes 
its physical value of 27.46(44) [10]. The preliminary value of is taken from Ref. [23]. Note 
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Figure 3: Linear extrapolations of 04, y/ioXM, ^ to the physical point (vertical red bars). The horizontal 
red bars indicate the corresponding experimental value including errors (or respectively the theoretical value 
including errors in the case of combinations involving fo), the blue horizontal bar represent the respective 
lattice value with errors obtained by extrapolating the data to the physical point. 


that is then determined by our definition of the chiral trajectory through TrM = const. We 
remark that this definition is preliminary since we would hope to improve over the errors on that 
ratio given in Ref. [10] once the continuum limit can be taken. 

We first look at the chiral extrapolation of 04 , see the upper right plot in Fig. 3. Remember that 
04 has been used to match our chiral trajectory along the symmetric line. We can now ask whether 
our matching condition at that point results in an agreement at the physical point, i.e. 04 obtains its 
physical value at the physical point which is actually what is desired. Since m\ + m\/2 is constant 
to leading order along that trajectory and in addition Iq is also only weakly dependent on the quark 
mass (see upper left plot in Fig. 3) we do not expect a strong quark mass dependence of 04 , in fact 
this is evident from the upper right plot of Fig. 3 and we find agreemenf of 04 af fhe physical poinf 
wifh ifs physical value wifhin errors. Up fo cufoff effecfs, fhis means fhaf we have indeed chosen 
fhe correcf value of TrM for our chiral frajecfory. However, since we do nol know much abouf fhe 
cutoff dependence af fhe presenf sfage, if is inferesfing fo compare how differenf quanfifies mafch 
af fhe physical poinf because differenf quanfifies could have very differenf cufoff dependence. Lef 
us fake a closer look af and ^. The chiral exfrapolafions are shown in fhe lower plofs of 

Fig. 3. Wifhin errors we find again agreemenf for fhose quanfifies af fhe physical poinf wifh fhe 
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corresponding physical values. Again we would conclude that we have chosen the correct value of 
TrM for our chiral trajectory, and, in addition, the cutoff effects seem to be small. Note, however, 
that without knowing the explicit cutoff dependence of the quantities under investigation it is not 
really possible to disentangle cutoff effects and effects stemming from a mistuning of our chiral 
trajectory, i.e. a wrong value of TrM. This will be clarified by a future study of our ensembles at 
different lattice spacings which will also enable an independent scale setting. 

We remark that the conclusions we have drawn above, of course, depend on the actual value 
and error on to which is a result from a lattice study [ 11 ] as to is not experimentally accessible. 
We therefore also plan to extract to from our simulations and once the continuum limit can be 
performed we expect to deliver an independent cross-check on that quantity. 

Finally, coming back to the issue of scale setting, we remark that all quantities discussed in this 
section can be used for scale setting. As we have found that all these quantities agree well at the 
physical point with their corresponding physical value it is evident that the lattice spacings extracted 
by the quantities investigated here will results in similar values. As our results are still preliminary 
we will refrain from quoting numbers but will return to this topic in a different publication. 

5. Summary and Outlook 

We have presented first results on the baryon octet spectrum from lattice simulations with 2+1 
flavors and open boundaries which has been generated within the CLS collaboration. Using leading 
order chiral perturbation theory we have chirally extrapolated results down to the physical point for 
N,L,E at pion masses ranging from ss 220 — 420 MeV at a lattice spacing a ^ 0.086 fm. We 
have found nice agreement with the literature. In addition, we have extrapolated %, (p 4 , y/toXj^i, 
■jf to the physical point and have checked that 04 , which is used to match our chiral trajectories, 
agrees at this point with its physical value. Furthermore, we have found that -y/to^Ov and ^ also 
agree with their corresponding physical values at this point suggesting that cutoff effects may be 
small and that the chiral trajectory has been chosen appropriate. A more detailed analysis including 
cutoff effects will follow in the future. 
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